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Abstract 

Motivated by the recent Kitaev's K-theory analysis of topological insulators and superconductors, 
we adopt the same framework to study the topological phase structure of Jackiw-Rebbi model 
in 3+1 dimensions. According to the K-theory analysis based on the properties of the charge 
conjugation and time reversal symmetries, we classify the topological phases of the model. In 
particular, we find that there exist Z Majorana zero-modes hosted by the hedgehogs/t'Hooft- 
Polyakov monopoles, if the model has a T 2 = 1 time reversal symmetry. Guided by the K-theory 
results, we then explicitly show that a single Majorana zero mode solution exists for the SU(2) 
doublet fermions in some co-dimensional one planes of the mass parameter space. It turns out 
we can see the existence of none or a single zero mode when the fermion doublet is only two. 
We then take a step further to consider four-fermion case and find there can be zero, one or two 
normalizable zero mode in some particular choices of mass matrices. Our results also indicate 
that a single normalizable Majorana zero mode can be compatible with the cancellation of SU(2) 
Witten anomaly. 
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I. INTRODUCTION 

In the recent years, Majorana zero-modes in condensed matter systems have attracted 
much attention due to its topological and/or symmetry protected degeneracy (similar to the 
topological ground state degeneracy of topological ordered states). The degeneracy leads 
to topological and/or symmetry protected non-Abelian geometric phases (up to an overall 
phase) as we braid the defects that carry Majorana zero-modes. Note that the Hilbert space 
of these zero-modes are highly entangled |T], implying the underlying topological nature 
against local perturbations (that may repect certain symmetry in the symmetry protected 
cases). Thus, these Majorana zero-modes, if exist, can help to realize the robust quantum 
computation through protected non-Abelian geometric phases. (We would like to point out 
that those protected non-Abelian geometric phases do not correspond to the standard non- 
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Abelian statistics. This is because the defects that carry Majorana zero-modes have a long 
range interaction and the overall phase of the non- Abelian geometric phases may depend the 
exchange paths. So in this paper, we will call such non-Abelian geometric phases, defined 
up to an overall phase, projective non-Abelian statistics.) . 

These fermionic zero-modes appear near the core of some topological defects, which pre- 
serve either charge conjugation (particle-hole) symmetry or time reversal symmetry. Espe- 
cially, such objects in 3+1 dimensions have been proposed by Teo and Kane in [2J. It then 
opens the possibility to realize the non-trivial projective non-Abelian statistics, the so called 
projective ribbon statistic in [3] in 3+1 dimensions. Note that, unlike the standard non- 
Abelian statistics, the non-trivial projective non-Abelian statistics can exist beyond 2+1 
dimensions. Moreover, the more general classifications of all possible Majorana zero-modes 
in the disordered systems of any dimensions has been done in [3H9] based on K-theory anal- 
ysis in the context of the topological insulators and superconductors [TUHT5] . The above 
mentioned 3+1 dimensional Majorana zero modes are one of the examples in this general 
classification schemes. 

On the other hand, in the community of high energy physics the relations between 
fermionic zero modes and topological defects have been well-studied, and it results in the 
connection with the index theorem [TB]- The pioneer works relevant for our discussions here 
are the Dirac zero modes around 3+1 dimensional monopole (or hedgehog) by Jackiw and 
Rebbi p2] and around 2+1 dimensional vortex by Jackiw and Rossi [18]. The zero modes 
of Majorana fermions, however, are only found in the cases of a kink in one dimension [T7j 
and a vortex background in two dimensions [18] . It is then interesting to see if there exist 
possible Majorana zero modes in the Jackiw-Rebbi models, which are similar to the ones in 
Jackiw-Rossi model, as long as the particle-hole symmetry or the time reversal symmetry 
are preserved by the background topological defects. Once this is possible, it is also inter- 
esting to explore the possible connection between the models in high energy physics and the 
models for topological insulators and superconductors. 

A general consideration along the aforementioned direction was first given by Hofava 
in |19j . trying to link the stability of Fermi liquids to the one of D-branes. Some works 
along this direction in a more specific context have been done in [2QH22] . In these works, the 
relation between Jackiw-Rebbi or Jackiw-Rossi models and the Bogoliubov-de Gennes (BdG) 
model for topological superconductor is explored and explicitly constructed. However, the 
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role of the SU(2) internal symmetry in the Jackiw-Rebbi model is not fully addressed in 
these works because the corresponding internal symmetry is not obviously implemented in 
the BdG model. In this paper, we will further explore this issue by studying the role of 
SU(2) symmetry in the charge conjugation and time reversal symmetry operations. In this 
way, we can find out the topological nature of the fermionic zero modes in the framework of 
K-theory analysis, and then explore the topological phase structures. 

In general, the topological defects such as monopoles can be realized and classified by the 
pattern of spontaneous symmetry breaking. For example, the monopole in the Jackiw-Rebbi 
model is the 't Hooft-Polyakov one, in which the SU(2) gauge symmetry is broken to U(l) 
by the Higgs vev of the hedgehog pattern. Thus, the Majorana zero modes, if exist, can 
be understood as the low energy phenomenon due to the background of topological defects. 
One would hope that the setting of the theory is also well-defined at high energy regime. 
This concern has been recently raised by McGreevy and Swingle in [23] for the two-flavored 
Jackiw-Rebbi model with nonzero Dirac masses. In this case, the theory of odd number of 
Weyl fermions will be plagued by the Witten anomaly |24j . However, as shown in [23] it 
is usually difficult to obtain odd number of normalizable Majorana zero modes out of even 
number Weyl fermions so that the expected projective non-Abelian statistics is in vain. One 
naive way to bypass this no-go is simply not gauging the internal SU(2) symmetry. Though 
the presence of the normalizable Majorana zero modes will not be affected by ungauging the 
internal symmetry, the relative motion of the hedgehog defects on which the Majorana zero 
modes live will require infinite energy. The reason is that without the cancellation from the 
gauge fields, the potential between the defects due to the Higgs scalars is confining. If so, 
these Majorana zero modes will be non-dynamical, not mentioning the nontrivial statistics. 

It seems that it is impossible to lift the above no-go result for the 3+1 dimensional Ma- 
jorana zero modes with nontrivial statistics. However, we will see that in the analysis of 
[23] the charge conjugation and time reversal symmetries are not imposed. In this paper we 
will re-examine the issue by requiring the charge conjugation and time reversal symmetries, 
which then restrict the flavor mass matrices. It turns out that these discrete symmetries are 
essential to have Majorana zero modes classified by integer Z according to the K-theory anal- 
ysis. Indeed, we find that the no-go theorem is lifted by the existence of 0, 1, 2 normalizable 
Majorana zero modes for some ranges of mass parameters in the 4-flavored Jackiw-Rebbi 
model, which is UV-complete without Witten anomaly. 
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In general, the symmetry of the point defects determine the possible structures of the 
zero-modes on the defect. In this paper, we consider the defects that have a time reversal 
symmetry (generated by T) and Z 2 symmetry (generated by C). The Z 2 symmetry can be 
realized as a charge conjugation symmetry. If T and C satisfy T 2 = C 2 = 1 and TC = CT, 
the symmetry group of the defect is given by G+ + (T,C) (using the notation of [?]). From 
the table VII in [7], we see that the different structures of the zero- modes are classified by 
Z. In this paper, we show that the different structures simply correspond to the number of 
Majorana zero-modes. 

If instead T and C satisfy T 2 = (—) Nf , C 2 = 1 and TC = CT (where N F is the total 
fermion number operator), then the symmetry group of the defect is given by Gl + (T, C). 
From the table VII in [7j, we find that the different structures of the zero-modes are classified 
by Z2, corresponding to even or odd numbers of Majorana zero-modes. 

Our paper will be organized as follows. In the next section we introduce the Jackiw-Rebbi 
model with many fermion species/flavors and their charge conjugation (C), time-reversal (T) 
and parity (P) symmetries. These symmetries also provide a starting point for the review of 
K-theory analysis on classifications of the topological insulators/superconductors in section 



III In section IV we solve the Majorana zero modes for the SU(2) doublet fermions directly 



from the equations of motion with two and four flavors. We find that there exists a single 
normalizable Majorana zero mode in both cases. This lifts the no-go result discussed in [23J. 
We also delineate the phase boundary in the mass parameter space. Finally we conclude 
with some comments in section |Vj Technical details about the general choices of the mass 
matrices's parametrization in four-flavor case will be given in Appendix. 



II. THE JACKIW-REBBI MODEL WITH MULTI-FLAVORED FERMIONS 

In this section we will introduce the models for classifying the topologically ordered phases 
in the context of high energy physics, and also later for finding the Majorana zero modes lo- 
calized on the point defect. These Majorana zero modes are argued to obey the non-Abelian 
statistics [31 [23]. The classification of the topologically ordered phases is relevant to finding 
the fermionic zero modes on the point defects, as suggested by the bulk/edge correspondence 
of the topological insulators/superconductors [15]. This correspondence relates the pattern 
of the edge gapless excitations of a topological insulator/superconductor to the one of the 
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bulk gapped excitations, for more conceptual discussion see [25J. 
A. The Jackiw-Rebbi model and the discrete symmetries 

As shown for examples in [5HHJ |25], the classifications of the topological insulators or 
superconductors are tied to the discrete symmetries of the fermionic systems such as the 
charge conjugation (C), parity (P) and time reversal (T) symmetry. To further discuss the 
role of discrete symmetries and be more specific, we first introduce the Jackiw-Rebbi model, 
which describes Dirac fermions coupled to a gauge field and a scalar field in the SU(2) 
fundamental/adjoint representation. The Lagrangian density is 

C = ^prfD^! - Auipi§ a T a ipj - mu^rtj , (2.1) 

where = — igA^T 0- is the covariant derivative with A ^ the gauge field, T a 's are 
the SU(2) generators and $ a, s are the scalar fields. The space-time indices are denoted by 
/i, v = 0, 1, 2, 3, the SU(2) indices by a, b = 1, 2, 3 and the flavor indices by J, J = 1, 2, • • • , N. 
The flavor matrices A and m are hermitian. Our convention for gamma matrices is the same 
as the one in [22] , and we will show how to transform this chiral representation to a real 
representation of gamma matrices later. In this study, we are interested mainly in classifying 
the bulk topologically ordered phases and then use the results as the guidance in searching 
for the normalizable fermionic zero modes on the point defect. Since we are only interested 
in the pattern of these zero modes not their relative dynamics, we will then ignore the gauge 
fields which play no role for our consideration. 

The Hamiltonian derived from the above Lagrangian (without gauge field) is 

H = ^hip= ilh (-ia k d k ® l nm ® 1 LJ + (5 ® $ a T n a m ® A u + P ® ® m 7J ) </w , 

= ^(-i)(a k d k + M)il; , (2.2) 

where /3 = 7 , a k = 7°7 fc for k = 1,2, 3. This form of the Hamiltonian can be understood 
as an effective Hamiltonian near the Dirac cone of some electron systems in condensed 
matters. The constant mass matrix M is used to gap the system to obtain topological 
insulators/superconductors provided that there are no zero eigenvalues for M. If the scalar 
$ a is topologically non-trivial, there may exist some normalizable zero mode solutions which 
satisfy: 

h(x)il> (x) = 0. (2.3) 
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Later, we will first classify the these zero modes by K-theory analysis and then use the 



results as the guidance to solve (2.3) for the zero modes explicitly. 

To be specific, in this paper we will consider a particular kind of scalar condensate, 
the hedgehog configurations as considered in the Jackiw-Rebbi model [T7J [23], for which 
$ a (r) = r a (j)(r )Q The core of this point-like defect may host various patterns of normalizable 
fermionic zero modes depending on what kind of the discrete symmetries or some U(l) 



symmetries are preserved by the Hamiltonian (2.2). 

Three discrete symmetries can be realized in the Jackiw-Rebbi model: parity P, charge 
conjugation C and time-reversal symmetry T ^ They are defined by (up to a phase) 

Pip{t, x)p- x = U p tp(t, -x), (2.4) 
CV(t,£)C'- 1 = Uj>*(t,x), (2.5) 
f^(t,x)f- 1 = U t i){-t,x). (2.6) 

The U c , U t , and U p are unitary matrices, which can be written in the chiral representation 
as following: 



for SU(2) doublet, i.e., T a = \r 



U c = H7V) ® (^ 2 ) ® 1, (2.7a) 
U t = (7V7 5 ) ® (^ 2 ) g) 1, (2.7b) 
U p = (7°) ® 1 ® 1, (2.7c) 



where r a 's are Pauli matrices; 

• for SU(2) triplet, i.e., (T a ) nm = ie nav ,. 



U c = H7 2 ) ® 1 ® 1, (2.8a) 
f/ t = ( 7 S 3 ) ® 1 ® 1, (2.8b) 
f/ p = ( 7 °) ® 1 ® 1. (2.8c) 



1 A typical profile for </>(r) is 0(r) = u tanh(fcor) with a constant parameter kg and the vacuum value of 

scalar field v. We will use this profile to perform the numerical calculations. 

2 Again, we assume these symmetry transformations do not involve the flavor space. 
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The invariance of the Hamiltonian (2.2) under these transformations give 



Ujh* 



-hXJi 



Ujh* = ml 



U p h = hU p , 



(2.9a) 
(2.9b) 
(2.9c) 



respectively. From (2.7), (2.8) and (2.9), the mass matrices A and m have to take the form: 



A* = A, m* 



m 



(2.10a) 



for SU(2) doublet, and 



A* = -A, m* = m 



(2.10b) 



for SU(2) triplet. On the other hand, as considering a real (Majorana) fermion, the Dirac 
spinor is subjected to the reality condition: 



tp c (t,x) = U c ip*(t,x) = ip(t,x). 



(2.11) 



The explicit form of U c in (2.7a) (also (2.8a)) will be needed later when solving the 



Majorana zero modes on the point defect, and the parity St = UtU% will be relevant for 
classifying Majorana zero modes on the point defect based on K-theory analysis. 



Some remarks are in order. Firstly, we shall emphasis that the specific form (2.2), rather 



than the general form commonly considered in [HI [7], is governed by the Lorentz invariance. 
The special structure of Dirac gamma matrices then assures that parity symmetry is given 



by (2.7c) (also (2.8c)) since the scalar is treated as a background which does not transform at 



all under symmetry operations. We also assume CTP theorem so that if U c is also specified, 
then Ut will be fixed accordingly. 

Secondly, the fermions considered here are fundamental particles, rather than the com- 
posite ones with non-trivial constituents or spin-orbital interaction which are commonly 
realized in condensed matter physics. The lack of nontrivial internal structure for ip then 
demands that U C U* = +1 if we would like to obtain a real fermion from a complex one with 



the supplementary condition (2.11), which indicates the importance of charge conjugation 
symmetry. This can be easily verified by observing (ip c ) c = U*U c ip = ip, which implies 



U C U* = +1. Our choices of U c in (2.7a) and (2.8a) are conformed to this condition. 



Thirdly, suppose we are in a real representation in which the condition (2.11) becomes 



■ifj = ip*. This implies that the Hamiltonian h is anti-symmetric and purely imaginary by 



Fermi statistics and hermiticity of h. This condition and (2.10) require that 



(2.12a) 



for SU(2) doublet, and 



(2.12b) 



for SU(2) triplet. We will show how to transform the chiral representation (2.7) into this 
"tilde" representation in the next subsection. 



B. The real representation gamma matrices and real fermions 



When discussing symmetry operators such as (2.7) and (2.8), we use the expression 
of chiral representation of Dirac gamma matrices. We did this in purpose since the chiral 
representation gives us a great simplification when it comes to finding the zero mode solutions 
[T7] [23] . It is, however, important to have a real representation (in the literal sense) for 
these symmetries if we adopt the K-theory to classify the bulk topologically ordered phases 
and the edge zero (gapless) modes for Jackiw-Rebbi model. In this subsection, we show 
how to transform the symmetry matrices from chiral (complex) representation into a real 
representation for SU(2) doublet case. This gives an explicit example for the more general 
discussion given in [25J. For SU(2) triplet fermions, the explicit construction has been given 
in [20]. 

For a SU(2) doublet fermion, the scalar field $ = $ a r a /2 is in a pseudo-real represen- 
tation, i.e., B^*B^ = — $ with the unitary matrix B satisfying BB* — — 1. The obvious 



choice is B = it 2 which explains the appearance of it 2 in (2.7). In order to impose a consis- 



tent reality condition on the spinor in (2.11), however, we also have to choose a pseudo-real 
representation for the Dirac gamma matrices in which there exists an unitary matrix A with 
A{^f^)*A^ = ±7 M and AA* = —1. We can then always write A as ia 2 ® 1.2x2 and combine A 
and B by a tensor product: 

U' c = ia 2 ®l 2x2 ®ir 2 . (2.13) 
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It is obvious that U' C U' C * = +1. For such a symmetric matrix U' c , we can always bring U' c to 
a new representation U c by a Takagi transformation in which U c = SU' C S T = lg x g [2S]- The 
only task is to find the unitary matrices S. 
Let us consider an unitary matrix S: 

' P_ -iP + 
iP+ -P_ 

P± is a projection operator with the properties 
P + + P_ = 1, P±P T = 0, Pj = P± and P± = P±. It is easy to check S is unitary and 
U c = SU' C S T indeed gives us an identity matrix, i.e., 



8x8 



(2.14a) 



where P± = |(1 ± aV). Note P| 



(2.14b) 



P_ -zP+ \ / (^ 2 )(zr 2 ) 
zP+ ~P ) \ (w 2 )(2T 2 

(P+) 2 + (P_) 2 

(P + ) 2 + (P_) 2 




•8xf 



Before we apply (2.14) to transform our chiral representation (2.7a) to the "combined" 



Majorana representation which mixes up the Lorentz and SU(2) gauge indices, we have to 



rotate the chiral representation of Dirac gamma matrices into the form (2.13). This can be 
easily done by 

1 1 



U' c = VU C V T \ with an unitary V = ^= 



® 1 



2x2- 



(2.15) 



-i i 



4x4 r 



where the subscript L indicates the matrix only act on Lorentz space. Hence we can trans- 
form our chiral representation to a Majorana representation by applying S and V sequentially 



u c = (sv)u c (svf, u t = (svyu t (sv)\ 

^=(5V)V, r = (SV)Y(SV)\ /i = 0,l,2,3. 



7° 



Note that 7^ and 7°r a , as promised in (2.12), are purely real and imaginary: 

-a 2 r 2 \ ^ _ / ia 3 r 2 

a 2 r 2 J' 7 ~ \ za 3 r 2 

ia 2 \ ^ -i^t 2 

to 2 / ' I -io 1 t 2 



7 2 



(2.16a) 
(2.16b) 

(2.17a) 
(2.17b) 
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-a 2 (r 2 r a 



t 2 (t")*) -i(r a + (r a )*) 
i(r a + (r a )*) 



CT 2/ T 2 T a_ T 2/ T ^*> 



(2.17c) 



Also, the reality condition (2.11) becomes ip c = ip* = ip 'm. this "tilde" representation. This 



can be verified as follows. We first write the 4-component Dirac spinor as 




(2.18a) 



By (2.11 ) we arrive 



?/, c = (_iy 7 5 ) (g) (ir 2 )^* = tp 



(ia 2 )(ir 2 )C 



(2.18b) 



On the other hand, from (2.16) we have 




(2.18c) 



where Re £ and Im £ are real and imaginary part of complex spinor £. From (2.18c), it is 
clear that ip* = if). This then completes our tasks in this subsection: seeking an unitary 
matrix (SV) to transform the gamma matrices to the real representation and the reality 
condition (2.11 ) becomes ip c = ip* = ip. However, since the equation of motion is much more 



simple and easy to deal with in the chiral representation, we will still work in this convention 



when solving the zero mode functions in Section IV 



III. CLASSIFICATION OF TOPOLOGICAL PHASES OF JACKIW-REBBI MODEL 

In this section we will classify the topological phases and the edge zero modes of the 
multi-flavored Jackiw-Rebbi model introduced above. This classification is based on the 
scheme of K-theory analysis first done by Kitaev in [6]. However, we will mainly follow the 
conventions given in [7j for our discussions, in which more general results in the context of 
condensed matters are considered. On the other hand, the more general results in context 
of high energy physics and AdS/CFT correspondence, see [25] . 
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The K-theory analysis is used to classify the gapped systems, i.e., deth 7^ 0, and the basic 
idea is what types of Clifford algebra can be formed for a given Hamiltonian h of massive 
free fermions in which all structures are fixed except the mass matrix M. The structure 
of these Clifford algebras yields a configuration space M. for M, whose homotopy groups 
7r n (A1)'s reflect the topological pattern of the bulk gapped/edge gapless excitations. 

To see the structure of the aforementioned structure of Clifford algebra, we first impose 
the "flat-band" condition h 2 = 1 since only the numbers of positive and negative eigenvalues 
but not their magnitudes are relevant as far as the topological property is concerned. This 
flat-band condition for Hamiltonian h then is translated to 

{a 1 ® 1, a j ® 1}= 25 ij ® 1, (3.1a) 
{a\M}=0, (3.1b) 
M 2 =-1®1. (3.1c) 

where the 1 here is the unit matrix containing gauge and flavor indices. Here we note that 



for a relativistic theory (3.1a) and (3.1b) are trivially satisfied owing to Lorentz invariance. 

Here an important point should be emphasized. When the flat-band condition h 2 = 1 is 
imposed, the space dependence of M (the hedgehog here) is turned off. In other words, we 
only focus on the matrix structure of M when it comes to the Kitaev's K-theory analysis. 
By doing this, the bulk/edge correspondence is implied manifestly, i.e., we are classifying 
the same mass matrix structure for bulk (no defect) and edge (defect present) Hamiltonian 
h. As we will see in the next section, the zero mode indeed lives on the defect while different 
structures of M gives us different number of zero modes. 

We first consider the classification of the topological phases for the complex fermions. For 



such cases, the flat-band condition (3.1) is all we need to classify the topologically ordered 



phases. This is because there is no way to form an enlarged Clifford algebra even if there 



exist some symmetries as we are going to discuss next. From (3.1), the configuration space 
Ai is determined by the complex Clifford algebra C7(3, 1) with all a*'s fixed, and is denoted 
by Ai = C3 in the convention of [7|. It then gives a trivial topological fermionic phase since 

MCs) = 0I3E5]. 

Now we turn to the case for the real fermions. Naively, if the discrete symmetries are 
not involved, then the configuration space is A4 = R® which again yields only trivial phase 
as 7r (R®) = 0. Here the pair (p,q) is referred to the signature of the Clifford algebra 
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Cl(p, q), and for the results of the homotopy groups of vr n (C p ) and 7r„(i2|), please check the 
summarized tables in [7J [23] . 

However, in this paper the model we consider does respect the charge conjugation and 



time reversal symmetries. In such a case, we apply these transformations (2.9) to (2.2) then 
yields the conditions on a"s and M as follows: 



a* U p = -U p a\ 


(3.2a) 


MU P = U P M, 


(3.2b) 


a l U c = U c (a l )\ 


(3.2c) 


MU C = U C M\ 


(3.2d) 


atyUt = -U t a\ 


(3.2e) 


M*U t = -U t M. 


(3.2f) 



where U p , U c and Ut are unitary matrices. 

To find the Clifford algebra for K-theory classification, some discussions are in order. Let 
us observe that it is impossible for U p to become one of the generators of an enlarged Clifford 
algebra since it is only anti-commutative with a 1 rather than both a 1 and M. Therefore, we 
will ignore parity hereafter. 



There are chances for either U c or Ut to join (3.1) to form a larger Clifford algebra. The 



only obstacle is the complex conjugate on ot and M as shown in (3.2). However, if there 



exists real representation for all the generators appear in (3.2), then the resultant Clifford 
algebra can be enlarged than the one without these symmetries. As shown in II B we indeed 
can find such a representation for Jackiw-Rebbi model. Moreover, in such representation 
U c — 1 so that it becomes trivial and plays no role in forming the enlarged Clifford algebra. 
However, we should caution the readers that there are more general consideration for which 
U c cannot be transformed into unity in the context of condensed matter physics. 

Moreover, in such a Majorana representation, all the gamma matrices are all either real 
or pure imaginary so that a l 's, M and U t constructed from them will also be all either real. 



In such a case, U t then anti-commutes with a l 's and M as seen from (3.2e) and (3.2f), it 
then yields an enlarged Clifford algebra. The structure of the resultant Clifford algebra will 
then only depend on the parity parameter 



s T 



u t ut 
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U 



(3.3) 



Fermions 


m 


A 


U e 




S T 


M 


bulk phases 


point defect 


Doublet 


Im 


Re 


(?7 2 7 5 ) <8> {it 2 ) 


(7 1 7 3 7 5 ) ® (it 2 ) 


+ 1 


Rl 





Z 


Triplet 


Re 


Im 


(-H 2 ) 


7 i 7 3 


-1 


R\ 


Z 


z 2 



TABLE I. The summary of symmetry operators and topological class for real fermions in Jackiw- 
Rebbi model. 

In the second equality we have assumed Ut to be real, so are a 1 's and M Q 

It is straightforward to show that st is invariant under Takagi transformation. Thus 
it is easy to check that st = 1 for SU(2) doublet, and st = — 1 for SU(2) triplet. For 
St = 1, in the convention of [7] the resultant Clifford algebra is C7(4, 1) so that M. = R\ 
and ^oiR-t) = yields only trivial bulk phase. On the other hand, for st = — 1 the resultant 
Clifford algebra is C/(3,2) so that Ai = R^. Using the result of Bott's periodicity theorem 
for K-theory, i.e., 



p p—q mod 8 > 



7r n (i?p) = 7T (-R° )., 



(3.4) 
(3.5) 



the bulk topological phases for real SU(2) triplet fermions are classified by vr (i?3) 
ir (R%) = Z. 



Moreover, we can also classify the edge Majorana zero modes with the help of (3.5). This 
is assured by the bulk/edge correspondence for the topological insulators/superconductors, 
which here is manifested as 



ir d -d b -i(M) 



(3.6) 



for the edge modes living on (4-dimensional defect. We are interested in the Majorana zero 
modes living on the point defect in this paper, c4 = and then the classification is given 
by ^(-R^) = Z for = 1 and ^(-Rg) = Z 2 for St = —1- These results are summarized in 
Table |T} Note that if there is no symmetry, the Majorana edge modes on the point defect is 
classified by ^2(^3) = Z2, however, in this paper we will not solve the zero modes for such 
a case. 

3 Due to the hcrmiticity of a 1 and anti-hermiticity of M, a l 's are symmetric and M is anti-symmetric if 
they are all real. 
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Before we close this section, we should note that the K-theory analysis discussed here can 
only guide us to find the number of zero modes hosted on the topological defect. However, 
it cannot tell us the detailed state space structure of these zero modes. This issue has been 
addressed in 2+1 Jackiw- Rossi model with the background consisting of integer number of 
vortices in [27] . A general formula for the dimensionality of state space is given explicitly and 
the mode operators with fermion number conservation are realized by a restricted Clifford 
algebra which shares some similarities with K-thoery approach. 



IV. MAJORANA ZERO MODES IN TWO- AND FOUR-FLAVORED JACKIW- 
REBBI MODEL 

In this section, we discuss the structure of mass matrices of Majorana fermions in the 
SU(2) doublet representation in Jackiw-Rebbi model and solve the coupled differential equa- 
tion numerically. We find that there exists either no normalizable zero mode or a single Ma- 
jorana zero mode in two-flavor case, which depends on the Dirac mass parameters. When 
the number of fermion species becomes four, we can find that the number of normalizable 
zero modes may be zero, one or two which also depends on the Dirac mass parameters. 



As we discussed in Section III, the number of normalizable zero modes hosted on the 
point defect is expected to be the same as the one indicated by K-theory analysis. Hence, 
the criteria to write down the appropriate mass matrices m and A in question are flat-band 



condition (3.1c) and charge conjugation symmetry (2.9a), which gives 



{A,m} = 0, (4.1a) 
A 2 + m 2 = l, (4.1b) 
A* = A, m* = -m. (4.1c) 



For the given scalar condensate $ a (r) (a hedgehog here) and proper boundary conditions 



for the mode function, (4.1) are all we need to find normalizable zero modes indicated by 



K-theory analysis in Section III 
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A. Equations of motion for the zero mode 



The equation of motion for the Majorana zero modes governed by the Hamiltonian (2.2) 



is 



(a k p k l nm l u + (3M nmIJ ) ip nJ 
with the reality condition on the Dirac field 

v> = r = u c r = 



(4.2a) 




(4.2b) 



as shown in (2.18b). 



To proceed, we notice that the SU(2) Lorentz x SU{2) interna i symmetry has been broken 
explicitly by the background scalar field to the diagonal SU(2), the subgroup which mixes 
spatial and internal symmetries [T7j. So we decompose the two component spinor £ as 

(6a) j = (gAb + gjVib)(ir 2 )ba (4.3) 
where i,j = 1,2 and a, b = 1,2 are SU(2) and spin indices respectively. With this decom- 



position (4.3), (4.2) becomes a pair coupled differential equations: 

d r hj + Ajj^-hj + SIjjKj = 



d r (r 2 h rI ) + Ajj^-I^Kj) - ftjj(r 2 hj) = 



where 



and define 



g(r) 



l/A h(r), g = fg r {r) 



^ 4 K 



r r 



ft 



mi 



(4.4a) 
(4.4b) 

(4.5) 
(4.6) 



so that Q is real for pure imaginary m. We made an ansatz g = fg r (r) because of the 



rotation symmetry of the background field. Note that (4.4) are the differential equations 
with real coefficients so that we can solve for the "Majorana" zero modes. 

In order to find a normal mode, we have to transform the coupled first order differential 
equations (4.4) into a set of decoupled second order differential equations^] First, we define: 



4 The derivation given here is closely related to [23j except there it imposes different condition on m and 
thus yields different pattern from ours for zero modes of two-flavor case. 
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Wi ^ {e r d sAHs)/2 )ij{hj{r) ^ 

W rI = (ef rdsA ^/ 2 )Ur 2 KAr)). 



Hereafter we ignore the flavor indices. Then we can rewrite (4.4b) as 

d r W r (r) = e^ dsM{s)/2 fl(r 2 h). 



Integrating (4.8a) gives 



(4.7a) 
(4.7b) 

(4.8a) 

(4.8b) 



In use of (4.7a) and (4.8b), we can rewrite (4.4a): 

r 2 e-^ dsM{s)/2 d r W(r)= -SI (r 2 h T ) 

= _o e -r«*«A*«/2 j r dfe^ dsA ^ 2 n(r 2 h(f)) 

= _ e +r ^(»)/2 n a J' d re-f dsA ^y 2 {f 2 h{f)) (4.9a) 

where in the third equality we have used the fact f2A + Afl = twice. Since we choose to 
work in the basis where both A and ft 2 are diagonal, differentiating (4.9a) arrives a set of 
decoupled second-order differential equations: 



(4.9b) 



where the Ajj and flu are diagonal elements of A and fl. In the following discussion, we 
will ignore the subscript i and use the unbold letters A and m to denote one of the diagonal 
elements for simplicity. 



The first thing we notice is that the solution of (4.9b) when the diagonal matrix element 
tt 2 = is 



h(r) 



e -fdsA4>(s)/2 | + 



~ +fdsA<t>(s) ( ± 



dre 



(4.10) 



Since we require the mode function is regular at the origin and the profile of 0(r) is <f>(r = 
0) = and <j)(r — > oo) = v (we set v — 1 throughout the paper for simplicity) FJ the 
regularity near the origin implies c\ = 0. Therefore, we can see that the solution h(r) is 
normalizable for A > as Q 2 = 0. 



A typical choice of <f>(r) is <f>(r) — v tanh^k^r) where kg is a constant. 
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For diagonal elements Q 2 7^ (f2? f = — < since f2 is a real anti-symmetric matrix), 



the solution for (4.9b) is 



c 



W( r -> 0) ~ - (e~ mr - e +mr ) (4.11) 

because we require W(r = 0) is regular and <ft(r = 0) = 0. For general r, the W(r) satisfies 
the differential equation 

r 2 W" + (2r - r 2 A0(r))W // - mVW = 0. (4.12) 
To derive the proper boundary condition for numerical calculation near r = 0, we first 



expand (4.11) near r ~ 0, i.e 

W(r « 1) ~ -2mc ( 1 + 
2 



2 2 
m r 



3! 



+ 



4 4 
m r 



5! 



W'(r << 1) ~ — m 3 cr 1 + 



m r 



10 



+ 



m r 
280" 



(4.13a) 



(4.13b) 



Then we use the polynomial expansion to solve (4.12) near r = and the solution should 
reproduce the result (4.13 ) when we turn off 4>(r). Here we use the profile of <p(r) = tanh(kor) 



where k Q is a constant, and the results are 

W(r) ~ a + air + a 2 r 2 + a 3 r 3 + a$r 4 + a$r 5 + 
W'(r) ~ ai + 2a 2 r + 3a 3 r 2 + 4a 4 r 3 + 5a 5 r 4 + ■ ■ 
0(r) ~ fc r 1 . 



Substitute (4.14) into (4.12) and compare the coefficients, we arrive 



_1 term: a\ = 0, 



Constant term: a 2 



a Q m 
6 



Ak + m 2 
r term: a 3 = — a% = 0, 

9 2A/c + ^ 2 m 2 (2Ak + m 2 
r term: a 4 = — a 2 = — 



20 

3 , 3A/c + m 2 

r term: as = as = 0. 

26 



120 



-a , 



(4.14a) 
(4.14b) 
(4.14c) 

(4.15a) 
(4.15b) 

(4.15c) 

(4.15d) 

(4.15e) 



So we take do = 1 and set c = ^ to match (4.13). Then the boundary conditions are 

W(r ->■ 0) - 



2 2 
m r 



m 2 (2Ak + m 2 ) 4 



6 



120 



m 2 r / 2Afc n + m 2 
H/'( r -> 0) — ( 1 + 



10 



(4.16a) 



(4.16b) 



Using the boundary condition (4.16), we can then adopt the numerical shooting method to 



find the normalizable zero mode solutions of the second-order differential equation (4.12). 
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B. Two-flavored case 



In two-flavor case, we can now parametrize our mass matrices to satisfy (4.1) as 



A = Aia 1 + A 2 a 3 , (4.17a) 
m = /iff 2 (4.17b) 

where Ai, A 2 and /i are real numbers with the supplement condition A 2 + A 2 . + /i 2 = 1. We 
make a rotation to the basis where A is diagonal, i.e., 



O 



( Ax A-A 2 \ 

V^ + CA-Aa) 2 yj A 2 +(A— A2) 2 
Ai — A2 — A 



(4.18a) 




A' = OAO T = I , (4.18b) 



m' = OmO T = m det(0) = j-} m (4.18c) 



where A = \J\\ + A 2 , > 0. In Fig. [T]we show the numerical result for the case /j ^ 0, which 
gives us no normalizable zero energy solution. If fi is tuned to be zero while keeping A the 



same, a single zero mode solution is popped out because of (4.10). The phase diagram is 
shown in Fig. [2j 

Here we should note a difference between the consideration in [23] and ours. We start 
with the hermitian mass matrices which are implied by the hermiticity of the Hamiltonian. 
On the other hand, in [23J they start from a generic two-Weyl-doublet model in which the 
hermiticity of their Hamiltonian is guaranteed by adding the hermitian conjugate rather 
than requiring the mass matrices to be hermitian. One may expect their mass matrices 
are more general than the ones considered here. It is, however, not true because they also 
require m to be hermitian after performing a Takagi transformation which makes A contain 
only real and positive diagonal entries. This actually imposes a constraint on the mass 
matrices and makes their model fall into different class from ours, so that there is no single 
normalizable zero mode found in their model. 
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FIG. 1. The profile of zero mode solutions of ( |4,4| ) for the mass matrices (4.18). The parameter 
choices are A = 0.6, /i = 0.4 and ko = 1. Left: h\(r). Right: h,2{r) 



FIG. 2. Phase diagram for (4.4) with mass matrix (4.18) of 2-flavored Majorana fermions. We 
assume A > 0. The normalizability of the zero modes only depends on the Dirac mass fx. The 
single normalizable zero mode appears only at the origin of this one dimensional phase space. 

C. Four- flavored case 

Now we would like to consider the case with four flavors of fermions, and we find that 
there exists a single normalizable zero mode in this case even for a Dirac mass matrix with 
some nonzero elements. We parametrize our mass matrices, which satisfy the requirement 
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(4.1c), by an outer product of two Pauli matrices: 

A =A 1 (a° g> a 1 ) + A 2 (a° <g> a 3 ) + A 3 (<7° g> a ) + A 4 (<x 2 ® a 2 ) 
+A 5 ((T 1 <g> <r°) + A 6 (a 1 <g> a 3 ) + A 7 ((x 1 <g> a 1 ) + A 8 (cx 3 ® a ) 
+A 9 ((T 3 <g> a 1 ) + Ai ((T 3 <g> (T 3 ), 

m =m 1 (o-° <g> a 2 ) + m 2 (o- 1 ® (j 2 ) + m 3 ((X 3 <g) a 2 ) 
+m 4 (o- 2 <g> a ) + m 5 (a 2 (g) a 1 ) + m 6 (a 2 (g> a 3 ). 



(4.19a) 
(4.19b) 



In addition, to meet the condition (4.1a) and choose a basis in which A is diagonal, (4.19) 
can then be further reduced to 



A =A 2 (o-° <g> a 3 ) + A 3 (o-° ® a ) + A 8 (o- 3 ® o-°) + Ai (a 3 ® a 2 
m =mi(o-° (g) a 2 ) + m 2 (cr 1 <8> a 2 ) + m 3 ((T 3 g) a 2 ) 
+m 4 ((T 2 (g) a ) + m 5 (cr 2 g) cr 1 ) + m 6 (cr 2 ® a 3 ). 

Note that the mass parameter space has been trimmed down largely. 



(4.20a) 
(4.20b) 



Next we consider one special case of mass matrices (Allb) which is one among three 



cases in (All), and leave the detailed discussions about how to derive these proper mass 
matrices in Appendix |A} The reason to choose this particular case for illustration is that the 
normalizability of the zero modes can be seen from a set of decoupled first-order differential 



equations rather than a set of second-order ones. For the other two cases in (All ), we have 
to solve the second order differential equations for the zero modes numerically. The special 
choices of mass matrices are 

A=A 2 (a°®o- 3 ) + A 8 ((T 3 g>o- ), (4.21a) 
m =m 2 (o- 1 <g> a 2 ) + m 5 (a 2 <g> a 1 ). (4.21b) 

In this basis, A and m aren diagonal and totally off-diagonal, respectively. After tedious 



manipulations, we are able to write down (4.4a) explicitly, i.e. 



d r h x + (A 2 + Ag) 


<P\r, 
2 




+ (m 2 + m 5 )h r 4 


= 0, 


(4.22a) 


<9 r /i 4 (A 2 + A 8 ) 


0(r) 
2 


-ri4 


- (m 2 + m 5 )h rl 


= 0, 


(4.22b) 


d r h 2 (A 2 A 8 ) 


(j){r) 
2 




- (m 2 - m 5 )h r3 


= o, 


(4.22c) 


d r h 3 + (A 2 - A 8 ) 


0(r) 
2 


h 


+ (m 2 - m 5 )h r2 


= 0. 


(4.22d) 
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FIG. 3. Phase diagram for (4.22) of 4-flavored Majorana fermions. We assume A2 > Ag > 0. 



The two normalizable zero mode solutions appear at the origin. Along two straight lines ni2 = 
±771,5, there exists only single normalizable solution. Other than these regions, there exists no 
normalizable solution. 



We can see that (4.22a), (4.22b) and (4.22c), (4.22d) are paired to form two copies of 



the two-flavor equations. For arbitrary nonzero values of | A2 1 7^ |Ag| and \m2\ 7^ |m.g | , there 



are no normalizable zero mode solutions in (4.22). However, there is a single zero mode 
solution when 1 777,2 1 = \™>§\ ( | A2 1 7^ | As | ) . For 7772 = 7775, the pair hi and do not have 
normalizable solutions while hs is normalizable and hi is non-normalizable if A2 — Ag > 
and vice versa. In the case of 7712 = —7775, the pair and /13 do not have normalizable 
solutions while h\ is normalizable and /7 4 is non-normalizable if A2 + Ag > and vice versa. 
Moreover, two normalizable solutions are present when m 2 = 7775 = 0. That is, we have 
found 0, 1, 2 normalizable zero modes depending the mass parameters. We show the phase 
diagram in Fig. [3j 

The phase diagram shows that this single Majorana zero mode exists only in a small range 
of mass parameter space - a co-dimensional one plane; in other words, a single Majorana 
zero-mode emerges only when some diagonal elements (but not all) of the square of Dirac 
mass m 2 equal to zero. The two normalizable solutions only appears at the origin of two- 
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dimensional (7712,7715) plane, which is similar to two-flavor case: a single zero mode solution 
appears only at the origin of one- dimensional space \x. This phase diagram can be viewed 
as a stack of two one-dimensional phase diagrams, each of which corresponds to the straight 



line 7772 = ±7775 respectively. This simply reflects the feature of (4.22), which is seen to be 
two copies of two-flavor equations. 

Moreover, the pattern of normalizable zero modes found above is in accordance with the 
result of K-theory analysis, namely, classified by ^(Rl) = Z. In contrast, in the 2-flavored 
case it is hard to tell if the pattern of the normalizable zero modes is classified by Z or Z2. 
This just reflects the fact that the K-theory classification needs large number of spectators 
to ensure stable results. 



V. CONCLUSIONS 

Conventionally Majorana zero-modes are realized by splitting a complex Dirac fermion 
into two real degrees of freedom, i.e., the Majorana zero modes are not fundamental and 
always come in pairs. An attempt to find a single real mode was made by McGreevy and 
Swingle in [23] , in which they try to realize a single normalizable zero mode emerging from 
two SU(2) doublet Weyl fermions to avoid Witten anomaly. They found that it is impossible 
to realize a single Majorana zero mode in an Witten anomaly-free theory. 

Despite their promising attempt, we turn to the K-theory analysis for help. By adopting 
flat-band condition and accommodating the discrete symmetries, which are essential ele- 



ments in K-thoery analysis, we constraint the mass matrix structure in (2.2). Furthermore, 
we introduce the Majorana fermions by imposing the reality condition on complex Dirac 
fields, which is commonly used in high energy physics. It turns out that the K-theory anal- 
ysis tells us that there could be any integer number of normalizable zero modes hosted by 
the hedgehog (Table I) and indeed we can see the existence of a single zero mode when the 
fermion doublet is only two. We then take a step further to consider four-fermion case and 
find there can be zero, one or two normalizable zero mode in some particular choices of mass 
matrices. 

In summary, a single Majorana zero mode can be realized in Jackiw-Rebbi model with 
some constraints given by the discrete symmetries. We believe there exists an index theorem 
to support this result since it is a result guided by K-theory. It will be an interesting future 
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project to find the index theorem behind. 
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Appendix A: Four-Flavored Mass Matrices 

With four flavors, we only need to extend our mass matrix to 4x4 matrices, we can also 
parametrize them by an outter product of two Pauli matrices: 



A =Ai(a° ® a 1 ) + A 2 (a° <g> a 3 ) + \ 3 (a° <g> a ) + \ 4 {a 2 ® a 2 ) 
+A 5 ((T 1 <g> <x°) + A 6 (cr 1 <g> a 3 ) + X^a 1 <g> a 1 ) + \ 8 (a 3 ® a ) 

+ A 9 ((T 3 (g)(T 1 ) + Aio((T 3 ®CT 3 ), 



(Ala) 



m =m 1 (a° <8> cx 2 ) + m 2 ((T 1 <8> o 2 ) + m 3 (a 3 <g> cr 2 ) 
+m A (a 2 g> a ) + m 5 ((T 2 <g> cr 1 ) + m 6 (a 2 <g> cr 3 ). 



(Alb) 



Applying constraint (4.1a) leads to six equations: 




A 3 m 3 + A 7 m 2 + A 8 mi - A 6 m 5 



Aim 5 + A 2 m 6 + A 3 m 4 + A 4 mi 



A 3 mi + A 4 m 4 + A 5 m 2 + A 8 m 3 



A 2 m 4 + A 3 m 6 + A 7 m 3 - A 9 m 2 



A 3 m 2 + A 5 mi - A 9 m 6 + Ai m 5 



Aim 4 + A 3 ??7.5 — \qiji2 + Aiom 2 



0. 



0. 







0. 



(A2d) 



(A2b) 



(A2a) 



(A2e) 



(A2c) 



(A2f) 
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Here we choose Ai = A4 = A5 = A6 = A7 = Ag = to make A diagonal and the set of 



equations (A2) becomes 

A 3 m 5 + Ai m 2 = 0, (A3a) 

A 3 m 2 + Ai m 5 = 0, (A3b) 

A 2 m 4 + A 3 m 6 = 0, (A3c) 

A 3 mi + A 8 m 3 = 0, (A3d) 

A 2 m6 + A 3 m 4 = 0, (A3e) 

A 3 m 3 + A 8 m! = 0. (A3f) 
The mass matrices now reduce to 

A =A 2 (a° ® a 3 ) + \ 3 (a° ® a ) + A 8 (a 3 ® a ) + \ 10 (a 3 ® a 3 ), (A4a) 
m =mi(o-° <S> cr 2 ) + m 2 (a 1 <g> a 2 ) + m 3 (cr 3 <g> a 2 ) 

+m 4 (cr 2 <g> a ) + m 5 (cr 2 <g> cr 1 ) + m 6 (cr 2 <8> cr 3 ). (A4b) 
Now the A is diagonal. We first notice the diagonal part of m 2 are 

m 2 ^ = (mi + m 3 ) 2 + (m 2 + m 5 ) 2 + (m 4 + m 6 ) 2 , (A5a) 

m 22 = ( m i + m 3) 2 + ( m 2 - m^) 2 + ( m 4 - ^6) 2 , (A5b) 

m 33 — ( m i ~ m 3) 2 + ( m 2 - m^) 2 + (w 4 + m 6 ) 2 , (A5c) 

m 2 4 = (mi - m 3 ) 2 + (m 2 + m 5 ) 2 + (m 4 - m 6 ) 2 . (A5d) 



The requirement of off-diagonal terms in (4.1b) to vanish gives us the constraints on the 
choice of m's: 

mim 2 = 0, (A6a) 

mim 4 = 0, (A6b) 

m 5 m 4 = 0, (A6c) 

m^ms = 0, (A6d) 

m 6 m 2 = 0, (A6e) 

m6m 3 = 0. (A6f) 



First we notice that we have to choose at least three of the m^'s to be zero to fulfill (A6): 

nil — m 5 — m 6 — 0, or (A7a) 
m 2 = m 3 = m 4 = 0. (A7b) 
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However, these two choices lead us to nowhere since all A's have to be zero in order to satisfy 



the condition (A3) 



Next we consider the case when four of the m^s are equal to zero. To satisfy (A6), only 
nine combinations are allowed: 



777,1 = Wl2 = ITI3 = ™4 


= 0, 


(A8a) 


777-1 — m 2 — m 3 — m 5 


= 0, 


(A8b) 


nil — m 2 — m 5 — m 6 


= 0, 


(A8c) 


mi = ni3 = m>4 = 


= 0, 


(A8d) 


nil — m 3 — m 5 — m 6 


= 0, 


(A8e) 


nil = rn^ = 7775 = 7776 


= 0, 


(A8f) 


m>2 = rn3 = ^14 = "^5 


= 0, 


(A8g) 


m 2 = m 3 = m 4 = m 6 


= 0, 


(A8h) 


7772 = "^4 = nls = 7776 


= 0. 


(A8i) 



Now let us stop the choosing process here for a while and back to ( A5 ) . A lesson from the 
two-flavor case is that the mode functions are both non-normalizable if m| 7^ 0. Hence 
finding a single normalizable zero mode relies on = and the signs of the eigenvalues 



of A. From (A5), we see that the mass parameters come in pairs: (7771,7773), (7772,7775) and 
(7774, me) and they contribute to diagonal terms as (mj + mj) 2 to two of them and (mj — rrij) 2 
to the other two. Therefore, we conclude that we can only have single normalizable zero 
mode when the m's vanish in pairs (with proper choice of A). This reduces the possible 



choices of m's (A8) to three: 



7771 = 7772 = "^3 = = 0, 

mi = m 3 = m 4 = m 6 = 0, 

7772 = "^4 = "7.5 = "^6 = 0. 



(A9a) 
(A9b) 
(A9c) 



Combing (A3) with (A9), we find the possible parameter choices for having only a single 



normalizable zero mode as follows: 



mi = m 2 = m 3 = m 5 = 0, A 2 = A 3 = 0, 

7771 = 777 3 = 777 4 = 7776 = 0, A 3 = Ai = 0, 

777 2 = 777.4 = 777 5 = 777-6 = 0, A 3 = A 8 = 0. 



(AlOa) 
(AlOb) 
(AlOc) 
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Besides these three choices, we only have zero or two normalizable zero mode solutions. The 

\ 



corresponding A and m matrices for (A9) are: 

/ Aio + A 8 



m= i 



A- 





A 8 - Aio 
-(A 8 + A 10 

y o o 

/ 



m 4 + m 6 

V m 4 — rriQ 

/ A 2 + A 8 
-(A 2 - A 8 ) 
A 2 - A 8 







-(A 8 -Aio)/ 
(m 4 + m 6 ) 

-(m 4 - m 6 ) 



\ 






m= i 



\ o 
/ o 




y m 2 + m 5 







-(A 2 + A 8 )/ 

— (m 2 + m 5 ) ' 

m 2 — Ms 

(m 2 — 777.5) 















/ 



m= % 



( A 2 + Aio 
A 2 + Aio 
A 2 - Aio 

V o 
/ o 

mi + 7713 

— (mi — 777.3) 

y mi — m 3 







-(A 2 - Aio) J 

-(mi + 777,3) 





(Alia) 



(Allb) 



(Allc) 
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